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highest  frequencies  and  the  >  30-MeV  proton  flux.  Comparisons  with  other 
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Using  Solar  Radio  Burst  Integrated  Fluxes  to 
Predict  Energetic  Proton  Flux  Increases 

I.  IVUIOPI  (TIO.N 

The  OiTurrence  ol  solar  flare  activity  is  accompanied  by  an  increase  in  the  flow 
'f  ener. mv  from  the  sir'  into  the  interplanetary  medium.  This  increased  energy  flow 
leave-  the  sun  m  the  form  of  various  emissions  such  as  radio  waves,  X-rays, 
energetic  fin -irons  and  protons,  and  alpha  particles.  The  energetic  protons,  if 
hev  reach  the  earth’s  environment,  make  their  presence  known  by  increased  radio 
wave  absorption  m  the  polar  ionosphere.  These  protons  also  disturb  some  human 
■ctivitii  s;  for  e.x  mure,  impeding  transpolar  radio  communications  | called  polar 
c  ip  ih.sorp+ion  1 1  •(  \>|,  lowering  tne  ionospheric  I J- region  height  wnicli  upsets  some 
radio  n  igatioi  t  unctions,  and  disrupting  both  manned  and  unmanned  space  ventures, 
lii  M  illie  ins*  dice-  ♦ne  sol  .1  ilari  induced  disturbances  are  of  sufficient  magnitude 
•a  disrupt  'dectria  i)  wit  distribution  svsiems  and  cause  problems  in  telephone  eir- 
tms  m  i  long  iisian  -e  pipeline  distribution  systems  such  as  the  Trans -Alaska 
mii'dme.  1  liese  an  1  oilier  disturbances  10  man's  activities  have  provided  an  ineen- 
1  c.  c  ui  si u  Iv  i  ne  sin,  ns  a.  liviiv,  and  the  impact  of  this  activity  on  the  terrestrial 

'  ic  a.  1; 1  .  >•  I  I  I  ill  |  .it  1 1  1:1  .!'l  A I  If  11  s*  1  >fig  i 

!.  1  ai  :  d  1*  ■  1 1 .  \\.  II.  at  ai  1  1 " 1 !  Sola  i1 -’1  e  r  rest  rial  Predictions  Proceedings,  \  id.  li, 

'A  irking  C  ci  mil  lien  iris  111  f  Reviews.  International  Solar -To  r  rest  rial  l’rc  - 
H  'lions  Th'Occc.irin's  ati  T  \\  oiTishcm.  Moulder,  (  omrado.  22-2 7  Anril  I"?!1. 


The  purpose  of  this  study  has  been  to  find  potential  relationships  between  tne 
maximum  energetic  proton  flux,  attributable  to  a  solar  flare,  and  tne  radio  emis¬ 
sions  from  that  flare.  Tne  effort  has  been  directed  to  finding  adjustments  to  the 
peak  particle  flux,  as  observed  by  earth  satellites,  for  the  position  of  the  flare  on 
tne  sun,  and  combinations  of  tne  radio  frequencies  observed,  in  order  to  obtain 
better  possible  relationships. 


2.  BACKGROUND 

In  relating  solar  radio  burst  phenomena  to  energetic  solar  particle  emissions, 

Castelli  et  al“  and  Castelli^  showed  that  a  particular  radio  burst  peak  flux  spectrum 

could  be  associated  with  the  occurrence  of  PCA's  of  2  dB  or  larger.  The  PCA's 

4  -7 

are  related  to  the  increase  in  the  flux  of  MeV  energetic  protons.  The  radio 
spectrum  was  found  to  have  a  U -shape  when  the  peak  radio  fluxes  in  the  meter  through 
centimeter  portions  of  the  radio  spectrum  were  plotted.  The  time  constraint  for 

g 

this  condition  was  that  all  peaks  must  occur  within  a  few  minutes  of  each  other. 

The  minimum  of  the  U-shaped  peak  flux  spectrum  occurs  in  the  decimeter  portion 

of  the  radio  burst  spectrum  with  the  peak  fluxes  rising  to  >  1000  sfu  (1  sfu  =  1  solar 

flux  unit  =  10  YV  m  ^  Hz  *)  in  the  meter  and  centimeter  portions  of  the  radio 

spectrum.  Figure  1  illustrates  two  U-shaped  solar  radio  burst  spectra  associated 

with  energetic  proton  increases  observed  in  the  vicinity  of  the  earth. 

Several  different  approaches  have  been  used  to  find  possible  wavs  to  relate  solar 

<1 

radio  hurst  parameters  and  energetic  particle  fluxes.  Straka"  studied  radio  burst 
oarameters  such  as  burst  peak  flux,  burst  mean  flux,  burst  duration  and  burst -time 
integrated  flux  or  energv  at  five  frequencies  in  the  600-  to  8800-MHz  frequency 
interval  and  related  them  to  the  peak  proton  flux  of  the  associated  particle  event. 

He  found  that  the  t ime -integrated  radio  fluxes  gave  the  best  correlations  with  the 
proton  fluxes.  He  also  found  tiiat  the  (jest  correlations  were  with  the  lower  fre¬ 
quencies  in  the  observation  interval  even  though  his  sample  sizes  mav  have  been  too 

<t 

small.  In  this  integration  of  the  radio  fluxes,  Straka  integrated  t  tit'  entire  flux 
increase  of  the  burst,  including  anv  assoc  iated  post  burst  increases.  Newell'** 
int  eg  rated  only  Mu*  st  a  rt  -to -maxi  mum  portion  of  Hie  radio  flux  me  la,  sc  at  H.c  same 
five  radio  frequencies.  He  found  tiiat  tne  best  correlations  between  integrated 
radio  fluxes  and  peak  proton  fluxes  were  at  t  e  higher  frequencies  <>t  U.e  trequein  \ 
band  studied,  lie  also  found  that  adjustment  of  the  observed  particle  flux  bo  im¬ 
position  of  toe  sou  fee  fin  re  on  tne  face  of  the  still  did  not  i  rep  i  ac.  ,  ♦  oe  ,  ■  -  r  r<  I  d  ion. 


(line  to  trie  large  number’  of  references  cited  ihme,  H;<  will  no*  b-  !r  b  d  ,, 
See  t{efe  ra  nees,  page  at.) 


I'liver-1  ‘  related  integrated  burst  fluxes  at  2.  8  GHz  to  PC  A  magnitudes;  ne  did  not 

find  tne  flare  position  correction  to  provide  any  improvement  in  tne  correlation. 

It  is  interesting  to  note  tnat  tne  2.  8-Gllz  frequency  used  in  the  Oliver11  study  falls 

in  the  decimeter  portion  of  tne  radio  spectrum;  this  is  the  minimum  of  the  U-shaped 

2 

burst  spectral  shape  found  by  Castelli  et  al  to  be  predictive  of  a  proton  event 

1 2 

having  occurred.  Akin'yan  et  al  studied  both  rising  portion  and  the  whole  burst 
integrated  radio  fluxes  at  frequencies  of  3  and  9  GHz.  In  their  study  they  also 
made  adjustments  of  the  proton  fluxes  for  the  heliocentric  position  of  the  flare. 
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Figure  1.  Peak  Radio  Flux 
Spectra  of  Two  Solar  Proton 
Flares,  Kach  Showing  the 
U-shaped  Spectral  Characteristic 
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(  liver,  K.  \V.  ( 1976)  Parent  - 1  lari'  emission  al  2,  3  (.11/  as  a  Pri-.lictor  of  the 
Peak  Absorption  of  Polar  Cap  Kvents.  N  F  I  .C  •'  TR  -201a,  Naval  Fleet  rnni .  s 
1  .aboratorv  Center ,  San  Diego,  l  al  i  forma  "2  132. 

tkin'van,  S.  T.  ,  Fomirhev,  V .  V .  ,  an  I  (  'lu  rlnk,  1.  M.  <1!'78l  Kstir  i'is  of  the 
mtensitv  of  solar  protons  from  the  in'i  cral  parameters  o!  r  i.  nnnvi  n  In 
bursts.  Geomagnel.  Ac  roil.  IMtNo.  -It: 


All  of  the  studies  mentioned  here  were  made  at  discrete  frequencies  in  the 
decimeter  and  centimeter  portions  of  the  radio  spectrum.  Groom  ^  also  considered 
several  different  time  and  flux  characteristics  of  solar  radio  bursts  such  as  peak 
radio  flux,  mean  radio  flux,  time  duration  of  the  burst,  and  mean  duration  of  the 
burst  as  predictors  of  the  peak  proton  fluxes.  He  found  the  mean  duration  of  the 
burst  to  be  the  best  predictor.  The  mean  duration  of  the  burst  is  the  time-integrated 
flux  or  energy  of  the  burst  normalized  to  its  peak  flux.  One  other  characteristic  of 
the  radio  burst  to  which  Groom  '  referred  was  the  total  energy  density,  E^,.  A 
representation  of  E,j.  may  be  found  by  means  of  a  double  integration,  first,  a  time- 
■n  tog  rat  ion  of  the  radio  flux  increases  is  performed  at  each  frequency  of  observation; 
men  the  ti me -integrated  fluxes  at  the  frequencies  of  observation  are  integrated 
across  the  observed  frequency  interval  to  give  the  total  energy  density,  E^.,  of  the 
radio  burst 

14 

Bakshi  and  Barron  studied  the  use  of  Er^.  as  a  predictor  of  >  10  MeV  proton 
flux  maxima  using  integrated  fluxes  at  frequencies  of  606,  1415,  2695,  4995  and 
8800  MHz.  Instead  of  time -integrating  the  entire  flux  increase  at  each  frequency, 
however,  they  time -integrated  only  the  portion  of  the  radio  burst  flux  increase  that 
contributed  to  the  U  -shape  spectrum  discussed  earlier.  In  some  instances  this 
t i me -integration  was  of  the  entire  burst  flux  increase;  in  others,  it  was  not.  fig¬ 
ures  2  and  8  illustrate  the  procedure  used.  Figure  2,  an  illustration  of  the  solar 
radio  burst  observed  at  the  Sagamore  Hill  Radio  Observatory  on  28  August  1966, 
shows  an  example  where  onlv  the  portion  of  the  burst  exhibiting  the  U-shaped  pro¬ 
file  in  the  maximum  power  spectrum  was  used  to  obtain  the  time -integrated  radio 
fluxes.  First,  the  time -integration  of  the  flux  increase  at  each  of  the  five  indicated 
frequencies  was  found.  At  each  frequency  the  trace  of  the  flux  increase  was  broken 
into  a  collection  (with  a  small  amount  of  smoothing)  of  plane  triangles,  rectangles, 
and  trapezoids  (Figure  2).  The  areas  of  these  plane  figures  were  then  found  and 
added  together  to  give  the  representation  of  the  time -integrated  flux  for  that  fre¬ 
quency.  The  five  time -integrated  fluxes  were  then  represented,  as  in  Figure  2,  to 
give  four  trapezoidal  figures.  The  areas  of  trapezoids  enclosed  by  the  adjacent 
frequencies  were  found,  and  finallv  the  four  areas  were  added  together  to  give  the 
E.j.  for  the  radio  burst  being  considered.  The  use  of  electronic  digitizing  circuits 
and  computers  should  make  it  possible  for  this  process  to  be  more  a<  <  uratelv  and 
rapidlv  performed  in  the  future.  These  U.K.'s  as  well  as  the  time -integrated  radio 
flux  at  the  five  frequencies,  were  then  correlated  with  the  peak  hourlv 
average  >  10  Me\  proton  flux  associated  with  the  radio  burst  being  considered. 

Id,  Groom,  I).  I..  (197!)  Forecasting  the  intensitv  of  solar  proton  events  from  the 
time  charaeteristies  of  solar  microwave  bursts.  Solar  Uhvs.  If (No.  1 ) :  1 7  1. 

- - - - -  V'W 

14.  Bakshi,  U.  ,  and  Barron,  \V  .  (  1 ' 1 7 B >  Fredirtion  of  the  Proton  Flux  Magnitudes 
from  Radio  Burst  Data,  \  I  G  R  1 . -  |  R -7  8 -0 100.  AD  AOaTOTl.  ' 
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of  Time -Integrated  Radio  Fluxes 
of  L'-Portion  of  211  August  1966 
Solar  Radio  Burst 
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Table  1  gives  sonic  of  the  correlations  found.  It  is  seen  that,  in  the  case  of  the 
integration  of  individual  frequencies,  the  highest  frequencies  gave  the  better  correla¬ 
tions.  With  the  exception  of  4995-XUIz  data,  the  correlations  of  the  individual  time- 
integrated  radio  fluxes  and  K  were  further  enhanced  by  the  adjustment  of  the  ob¬ 
served  [jeak  hourly  averaged  particle  fluxes  for  the  location  of  the  source  flare  on 

•j  \ 

the  lace  of  the  sun.  1  he  adjustment  factor  used  was  e  ‘  ,  where  A  is  the  longitudinal 
distance,  in  radians,  from  the  location  of  the  flare  to  the  solar  foot  point  of  the  mag¬ 
netic  field  lines  connecting  the  earth  to  the  sun  at  ~57°\V  solar  longitude. 

Barron  and  Bakslii  ’’  found  that  a  better  correlation  between  K.  and  the 

■>  a  2  A 

■  10  -Me\  peak  particle  flux  was  obtained  bv  using  e"‘  instead  of  e'  ‘  .  It  has  also 
2  A 

been  found  that  the  e  adjustment  further  improves  the  correlations  at  the  top 
three  individual  frequencies  in  the  frequency  interval  considered,  as  shown  in  the 
last  i  n bin  n  iif  Table  I . 

It  will  be  the  purpose  of  ibis  report  to  discuss  the  effect  of  adjusting  the 
lootpoint  position  for  different  solar  wind  speeds,  and  different  combinations  and'or 
weightings  of  the  time -integrated  radio  fluxes. 


Barron,  U.  R.  ,  and  Bakslii,  P.  (  1 ' ‘7 f 1 )  Application  of  Integrated  Radio  Burst 
Fluxes  to  the  Prediction  of  Solar  Fnergetic  Proton  Flux  Increases,  Solar- 
i  crrest  rial  Predictions  Proceedings,  \ol.  Ill:  Solar  Activity  Predict  ion’s, 
international  Solar  Terrestrial  Prediciinn  proceedings  and  Workshop, 
Boulder,  ('otorado,  26-27  April  197*'. 


national  Astronomical  Union,  and  various  observatory  reports. 


:i.l  Radio  Data 

The  radio  data  are  primarily  data  taken  at  frequencies  of  606.  1415,  '2695, 
4095,  and  8800  MHz  at  the  Sagamore  Hill  Radio  Observatory,  the  Manila  Observa¬ 
tory  in  the  Philippines,  and  the  Athens  Observatory  in  Greece.  These  observations 
provided  a  consistent  long-time  set  of  data  at  the  five  frequencies  mentioned.  The 

peak  I'lux  radio  burst  spectrum  of  all  of  the  events  considered  had  the  U-shaped 

•) 

spectrum  described  bv  Oastelli  et  al.  "  In  two  of  the  events,  1  and  20,  radio  data 
information  from  observatories  other  than  the  three  mentioned  were  used.  In 
event  number  1,  24  March  1066,  no  observations  were  made  at  the  live  indicated 
frequencies  at  the  time  the  flare  occurred,  ~  0225  I  T.  In  their  place  the  re¬ 
ported  values  at  500  MHz  from  (Iiraiso  Observatory  and  1000-.  2000-,  .1750-, 
and  9400-MHz  values  from  Nagova  Observatory  in  Japan  were  used.  Kvent  20 
was  viewed  bv  the  -V  hens  Observatory  whi  eh  does  not  observe  at  the  606  -  MIlZ  frequent' 
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TabU-  2.  Radio  ami  Proton  Data,  Radio  and  Solar  Wind  Data  (Contd) 


Radio  and  Solar 

Wind  Data 

Integration  of  T  (x*)'s  Across  the 
index  606-8800-MHz  l-rcquenev  Internal1 

V  c* 

solar  wind 

404/^  solar  wind 

1 

10 

6 

458 

(1) 

88 

J 

1  10 

6 

442 

(1) 

01 

;; 

17  a 

!» 

4  14 

(1) 

0. 

0  8 

4 

17  10 

0 

2  64 

(1) 

1. 

1  1 

12  1 

8 

640 
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0. 

63 

?, 

:s  i 

5 

.27  2 
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1. 

00 

V 

7  2 

2 
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7  1 

i\ 

t!4 
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0. 

86 

f»7 

7 
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59 
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01 
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454 

(1) 

0. 

80 

1 1 

1 1 1 
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0. 

88 
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(1) 

0. 

88 
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(1) 

0. 

M0 

!  1 

7 

188 

(  1  ) 

0. 

82 
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(1) 

0. 
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50 
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lieobh.vsieal  Data 

In  its  place,  thi-  radio  information  at  536  MHz  observed  at  ( indro  jov,  Czechoslovakia 
was  used.  The  inlcm  ated  radio  fluxes  lor  all  the  frequent  ies  of  event  1  and  the 
.536-MIiz  frequent  v  •>:  event  20  were  tabulated  from  the  event  mean  flux  multiplied 
hv  >he  event  duration  beeause  the  analog  irai  cs  of  these  events  were  not  available 
‘or  mteyralion  of  the  I  vpo  Oeserdied  earlier.  1  he  last  three  flares  m  the  table,  all 
>e.  urned  alter  the  1 07  fj  m  i  ni  mum  between  solar  eveles  20  and  2  1.  Twooftiiose 
dares,  20  and  21,  were  observed  ’()  oeeur  at  21°  N  solar  latitude,  with  ll  would 
m  ike  ilvm  probable  ohenomena  of  evele  2  1.  1  dare  lo  was  observed  at  07’  \  solar 


latitude,  meaninu  it  was  probablv  a  vestigial  event  of  the  20tll  sunspot  evele. 


3.2  Proton  Data 


The  particle  data  consists  primarily  of  the  peak  hourly  average  proton  fluxes 
at  '  10  MeV  and  >  30  MeV  as  viewed  by  various  Explorer  satellites.  The  first 
three  events  and  the  last  three  are  reconstructed  from  other  energy  channels,  as 
>  10  Me\  and  >  30  MeV  channels  were  not  available  when  these  events  occurred. 
The  particle  flux  values  used  were  adjusted  by  subtraction  of  background  particle 
flux  levels  and.  when  the  increase  under  consideration  occurred  on  top  of  a  prior 
increase,  any  flux  due  to  the  prior  increase  was  estimated  and  also  removed.  One 
additional  possible  limitation  on  the  observed  proton  flux  values  was  pointed  out 
tiv  t'room*'  where  lie  indicated  that  the  true  particle  flux  peak  may  not  be  well 
represented  bv  the  hourly  average  values  being  reported.  When  the  averaging 
process  is  considered,  it  will  be  realized  that  the  only  instance  where  the  average 
>f  a  group  of  values  will  equal  the  peak  or  highest  value  in  the  group  is  when  all 
values  in  the  group  are  equal  to  that  peak  or  highest  value. 


V.  (MRF.l.ATIONS 

Observations  of  solar  flare  phenomena  by  such  means  as  multifrequency  radio 
burst  measurements,  II-alpha  and  optical  observations,  and  X-ray  measurements 
all  provide  information  about  what  is  happening  in  the  vicinity  of  the  location  of  the 
flare.  The  radio  flux  increase  measurements  give  an  indication  of  protons  and 
electrons  that  have  been  accelerated  at  the  site  of  the  flare.  If  this  location  is  a 
distance  from  the  foot  point  of  the  magnetic  field  lines  connecting  the  sun  and  the 
earth,  there  is  going  to  lie  a  decrease  in  the  number  of  flare  accelerated  particles 
that  will  reach  the  sun-earth  magnetic  footpoint.  There  will,  therefore,  be  fewer 
particles  to  escape  into  the  interplanetary  medium  along  the  field  lines  to  the  vicinity 
>f  the  Earth.  Attempts  have  been  made  to  compensate  for  this  particle  loss  in 
various  prediction  studies,  as  mentioned  earlier. 

\  an  llollebeke  et  al*(’  noted  that  the  peak  dO-iMeY  particle  flux  observed  by 
earth  satellites  appeared  to  decrease  by  two  orders  of  magnitude  for  every  radian 
oetween  the  flare  site  and  what  thev  referred  to  as  the  "preferred-connection" 
region.  This  corresponds  to  the  "optimum  longitudinal  interval"  (01,1)  of 
Akin'van  et  al*1  but  is  not  necessarily  the  "fast  propagation  region"  of  Keinhard 

13.  \  an  llollebeke,  M.A.I..  Ma  Sung,  1..S.  ,  and  McDonald,  K.  15.  (  lf)75)  The  varia - 

lion  of  solar  proton  energy  spectra  and  size  distribution  with  helinlongitude. 
Solar  I’livs.  I  l(No.  1):  1H:>. 

-  w 

IV.  Akin'van,  S.  I.,  Eomichev,  Y.V.,  and  Chertok,  I.M.  (  1  !*77  )  Ueterminal  ion  of 
the  parameters  of  solar  protons  in  the  neighborhood  of  the  earth  from  radio 

bursts:  1.  Intensity  i  unction,  (leomagnet.  Aeron.  17(No.  D:ii. 
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1  H 

and  VVibberenz.  With  the  use  of  data  from  Table  2,  a  comparison  was  made  with 
the  Van  Hollebeke  et  al16  102A  coronal  gradient  attenuation  factor  to  see  if  some  other 
exponent  coefficient  would  provide  a  better  correlation  between  the  observed 
>  10-MeV  particle  flux  and  the  double -integrated  radio  flux.  The  adjustment  mul- 

n  A 

tiplier  used  with  the  peak  >  10-MeV  proton  flux  was  10  .  where  0  ‘  n  ?  2.  5  and  A 

was  the  solar  Longitudinal  angular  distance,  in  radians,  from  the  location  of  the 
flare  to  the  footpoint  of  the  magnetic  connection  between  the  sun  and  the  Earth, 
taken  to  be  at  57  °  W  solar  longitude.  This  is  approximately  the  center  of  the 
"preferred-connection"  region  located  in  the  20°  to  80°  W  solar  longitudinal  interval. 
With  A  determined  for  each  event,  the  problem  was  to  find  what  value  of  the 
coefficient  n  would  give  the  highest  correlation  between  the  adjusted  >  10-IVleV  proton 
neak  flux  and  frequency-time -integrated  radio  flux,  E^.. 

The  variates  correlated  were: 

X  =  Log  E-p,  the  independent  variable,  and 
n  \ 

Y  =  Log  (10  dependent  variable. 

I  is  the  peak  hourly  average  >  10-MeV  proton  flux. 

Figure  4  shows  the  variation  of  the  correlation  of  Y  with  X  for  different  values  of 
the  coefficient  n.  The  highest  value  of  r,  the  correlation  coefficient,  is  seen  to 
fall  in  the  interval  1.0  <  n  <'  1.5.  Solution  for  the  maximum  of  the  best  fit  parabola 
to  the  six  (n,  r)  pairs  gives  n  -  1.  274  and  r  =  0.  850.  With  use  of  this  value  of  n,  the 
adjusted  1  ^  fluxes  were  calculated  and  the  coefficients  of  the  best  fit  straight  line 
determined  to  give 


Y  =  1.  845  X  -  1.  107 


where  X  and  Y  are  as  given  here.  This  mav  also  be  written  as 


where  1  ■■  peak  satellite  observed  .»  1()-Me\  proton  flux.  Figure  5  illustrates  the 

poinls  from  Fable  2;  the  straight  line  is  given  bv  the  expression  above. 

The  time-integrated  radio  fluxes  at  each  of  the  five  frequencies  being  considered 
m  ibis  stui Iv  were  also  correlated  with  the  >  10-MeV  proton  flux  adjusted  as  men- 
t  Mined.  The  results  of  these  correlations  and  t he  coefficients  of  the  best -fit  straight- 
line  determined  from  the  data  set  in  'Fable  2  are  given  in  Fable  2,  along  with  the 


111.  lb  iriliar  i,  k.  ,  and  Wibberenz,  (I.  (1074)  Propagation  of  flare  protons  in  the 
solar  atmosphe  re.  Solar  I’hvs.  26(\o.  2):472. 
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aforementioned  K,j.  correlations.  It  is  noted  that  the  correlations  are  higher  at 
the  higher  three  frequencies  than  at  the  lower-  two.  In  addition,  the  combination 
of  the  time -integrated  radio  fluxes.  H^.,  correlates  better  with  the  adjusted  proton 
fluxes  than  do  any  of  the  individual  frequencies. 
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1.1  Solar  Wiml  Speed  Corn'l  l  ion 

The  discussion  thus  far  lias  assumed  that  the  foot  point  of  the  magnetic  field 
lines  connecting  the  sun  and  the  Earth  was  fixed  at  approximately  57  degrees 
(1  radian)  west  of  the  solar  central  meridian.  In  reality,  the  footpoint  of  the  eon¬ 
neeting  field  lines  will  varv  in  position  relative  to  the  solar  central  meridian  with 

19 

changes  in  the  speed  ot  the  solar  wind.  The  higher  the  velocity  of  the  solar  wind, 
the  closer  to  the  solar  central  meridian  will  be  the  footpoint  of  the  magnetic  field 
lines  which  are  frozen  into  the  solar  wind  plasma.  The  relationship  between  the 
solar  wind  speed  V  in  kilometers  per  second  observed  at  the  Earth  and  the 

longitudinal  position.  0.  ,  in  radians,  west  of  the  solar  central  meridian  as  given 
•>0  •' 

by  Smart  and  Shea“  is; 


404 


This  expression  was  used  to  find  the  different  footpoint  positions  for  the  events 
listed  in  Table  2.  The  solar  wind  speeds  used  were  the  closest  speeds  preceding 

2  l 

the  time  of  the  solar  flare  involved,  as  given  bv  King.  The  pre flare  solar  wind 
speed  was  used  because  present  considerations  are  for  the  purpose  of  predicting 
the  proton  flux  increase  from  a  solar  flare.  If  the  prediction  were  held  up  for  the 
changed  \'  ...  due  to  the  flare,  the  energetic  particles  mav  already  have  arrived, 
making  the  prediction  unnecessary.  ““ 

The  adjustment  factor  used  in  this  case  is  enl\  where  n  is  a  constant  to  be 
determined  and  B  is  the  angular  distance,  in  radians,  between  the  longitude  of  the 
flare  and  longitude  of  the  footpoint  found  as  described  earlier.  The  procedure  was 
the  same  as  in  the  previous  section:  to  find  the  value  of  n  which  produced  the  best 
correlation,  and  then  the  coefficients  of  the  best  fit  straight  line  to  the  data  set. 
Values  of  n  were  obtained  from  correlations  of  the  total  energy  density,  E  ,  in  the 
606  to  8800-MHz  frequency  interval  with  adjusted  proton  fluxes  at  both  lO-.MoY 
and  ■  30-MeV  energy  levels.  Figures  6a  and  6b  snow  the  correlation  coefficients 
obtained  for  different  values  of  n  for  the  >  10-AIeY  and  20-MeY  protons,  re¬ 
spectively.  The  values  of  n  for  maximum  correlation  from  the  two  plots  are; 

njg  =  2,  682  for  >  10-MeY  protons. 


and 


n,^Q  2.294  for  >  20-MeY  protons. 


(Due  to  the  number  of  References  cited  above,  tliev  will  not  hi'  listo.l  here. 
See  References,  page  2  1.) 
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Figure  6.  Variation  of  the 
Correlation  of  Log  ( e n i:iI )  With 

Log  Eq-  for  Different  Values  of  the 
Coefficient  n.  B  is  the  distance, 
in  radians,  from  the  flare  site  to 
the  solar  footpoint  of  the  magnetic 
Earth-Sun  connection  as  determined 
from  the  solar  wind  velocity. 

Figure  6a:  K  =  10,  that  is, 

>  10-MeV  protons 
Figure  6b:  K  =  30,  that  is 

>  30-MeV  protons.  The  points  of 
maximum  correlation  are  indicated 


By  incorporation  of  these  n-values  in  the  e11  *  adjustment  factor,  straight-line  fits 
were  made  to  the  adjusted  proton  fluxes  and  the  E.^'s  of  the  burst  and  the  five  time- 
integrated  radio  fluxes.  Table  4  gives  the  results  of  these  straight-line  fits.  The 
three  higher  frequencies  are  again  seen  to  correlate  be“er  than  the  two  lower  fre¬ 
quencies.  It  is  also  interesting  to  note  that  the  >  30  MeV  correlations  are  slightly 
better  at  the  top  three  frequencies  than  the  >  10  MeV'  correlations.  Plots  of  the 
best  fit  straight  lines  and  the  points  used  to  determine  them  are  given  in  Figure  7 
and  f  igure  fl  for  ■  10  MeV  protons  and  >  30  MeV  protons,  respectively.  The 
equations  of  the  two  straight  lines  are 


KM 


>  10  MeV 

V  =  1.  696  X  -  0.  991,  where 

X  =  Log  Et  (606-8600  MHz)  and 

..  ,  2.  682  H 

\  =  Log  1 , Qe 


>  30  MeV 

Y  =  1.  642  X  -  1.413,  where 
X  =  Log  E,„  (606-8800  MHz)  and 

x 

v  7  ,  3.  294  B 

Y  =  Log  I  ne 


where  B  is  the  solar  longitudinal  angular  distance,  in  radians,  from  the  location  of 
the  flare  involved  to  the  magnetic  footpoint  location  as  determined  from  404/V™.. 

An  observation  on  the  results  obtained  in  this  consideration  is  that  the  value 
n  =  3  used  by  Smart  and  Shea,  is  between  the  values  of  n  found  in  this  investiga¬ 
tion  for  the  >  10  Met  proton  fluxes  and  >  30  MeV  proton  fluxes;  for  this  reason 
n  =  3  can  be  used  with  little  degradation  of  results  for  the  correction  factor  in  both 
energy  channels.  It  is  also  a  simpler  number  to  use  in  anv  computation  scheme. 


Table  4.  Kadio  and  Proton  Correlations  (magnetic  footpoints  determined 
from  solar  wind  speeds) 


frequency 

(Mlz) 

Correlation  Straight-line  (Y  = 

Coefficient  r  b(intercept) 

mX  +  b)  Coefficients 
m(slope) 

>  10  MeV 

protons 

8800 

0.  802 

0.  4  90 

1.516 

4  995 

0.  8  1 5 

0.  529 

1.  742 

2695 

0.  780 

1.  125 

1.  64  9 

1 4  1  5 

0.  76  1 

1.  688 

1.  266 

606 

0.  7  03 

1.781 

0.  823 

606-8800 

0.  86  1 

-0. 091 

1 .  606 

'•  30  Me\ 

protons 

8.800 

0.  825 

-0.  064 

1.  5  29 

1995 

0.  826 

0.  007 

1.  7  30 

26  95 

0.  7  83 

0.  6  10 

1. 6  25 

14  15 

0.  ,  18 

1.  280 

1.  140 

6  06 

0.  646 

1.  3  20 

0.  74  2 
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0.  850 

-1.4  13 
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Figure  7.  Peak  Proton  Fluxes,  Multiplied  by 

\ppropriate  en**.  Plotted  Against  E^,.  The 
straight  line  best  fit  to  the  points  is1  shown. 

>  10  MeV  protons 

1.2  Multivariate  Prediction  Study 

The  results  of  the  correlations  performed  in  the  two  preceding  sections  indi¬ 
cated  that  the  best  relationship  between  integrated  radio  Ruxes  and  adjusted  ener¬ 
getic  proton  fluxes  were  obtained  when  the  total  energy  density  E^.  in  the  observed 
frequency  band  was  used, 

Results  also  indicated  that  in  taking  the  five  frequencies  involved  individually 
the  upper  three  had  better  correlations  with  the  proton  fluxes,  in  this  section, 
various  multivariate  combinations  of  the  integrated  fluxes  at  these  frequencies  will 
be  made  to  provide  linear  fits  to  the  proton  data. 


-  30  MeV  PROTONS 


10  >00  1000  10,000 

i  ■  (606-8800  MHz) 


Figure  8.  Peak  Proton  Fluxes,  Multiplied  by 

Appropriate  e11^,  Plotted  Against  K^..  The 
straight  line  best  lit  to  the  points  is  shown. 

>  8  0  MeV  protons 

4.8  Two  variiilc  Linear  Regression 

Two-variate  linear  regression  in  the  i'amiliar  straight-line  fit  for  finding  the 
coefficients  of  the  expression  V  =  m.\  *  b.  This  two-variate  linear  regression  was 
used  in  the  previous  sections  employing  onlv  time -integrated  radio  fluxes  at  dis¬ 
crete  frequencies  or  the  total  energy  densitv,  K.|, ,  as  tlie  independent  variable  and 
the  corrected  peak  proton  flux  as  the  dependent  variable.  Mere  we  will  investigate 
some  different  combinations  of  the  time -integrated  radio  fluxes  to  provide  various 
independent  variables. 

In  reference  to  f  igure  8,  'here  are  four  trapezoidal  areas  indicated,  two  of 
which  involve  the  lower  three  frequencies,  and  two  of  which  involve  Ihe  upper  tbrei 
frequencies  of  the  five  used  in  this  study.  The  frequency  261*5  MHz  is  seen  to  be 
included  iti  both  the  upper  and  lower  frequency  sets  of  t  rapezoids.  From  Figure  d, 


the  areas  of  trapezoids  1  and  2  were  combined  to  give  an  E^, -lower  and  the  areas 
of  trapezoids  3  and  4  were  combined  to  give  an  E^ -upper.  The  E^, -lower  and 
E^. -upper  for  all  the  events  in  Table  2  were  then  used  as  the  independent  variables 
to  obtain  the  coefficients  of  the  best  fit  linear  regressed  expression.  Again  the  fits 
were  made  to  the  logarithms  of  both  the  independent  and  dependent  variables.  Some 
fits  also  were  made  using  the  actual  values  of  the  data.  The  results  were  not  very 
different  as  far  as  the  correlations  between  the  variates  were  concerned,  but  the 
.  oet'ficients  obtained  were  much  larger  numbers  which  would  be  more  cumbersome 
u>  utilize  in  a  computing  scheme.  Table  5(a)  gives  the  results  of  the  straight  line 
fits  to  K.j. -lower  and  H.,.-upper  with  the  peak  adjusted  proton  fluxes  at  >  10  MeV 
and  >  30  MeV  energies. 


Table  5.  Radio  and  Proton  Correlations  (various  combinations 
of  time -integrated  radio  fluxes) 


(a)  Two-Variate  Linear  Regression,  Y  =  mX  +  b 


Low  Frequency 


High  Frequency 


>  10  MeV 


>  30  MeV 


>  10  MeV 


>30  MeV 


0.  7  52 
0.  57  1 


0.  826 
-0.  874 


-1.  397 


(b)  Three -Variate  Linear  Regression,  Y  a^  +  a^  Xj  *  XV 


Low  Frequency 


High  Frequency 


>  10  MeV 


30  MeV 


10  MeV 


>  30  MeV 


0.  792 


-0.  107 


-0. 446 


0.  842 


l.  1  Him-  Vilriiltc  Linear  Regression 

Instead  of  .  ombining  the  tw  >  areas,  1  and  2  tor  F^. -lower  and  3  and  4  for 
K.|. -upper,  in  Figure  3,  a  three -variate  linear  regression  lit  of  the  form 
V  n  •  ;ij  Xj  ■  \,t  vtns  performed.  The  variates  in  the  expression  .are 


I.oq  (adiusted  proton  flux) 

i  l.  )”  (Trapezoid  1  arrn  from  !•’  . -lower) 
or 

l  Lot;  (  I  rapozoid  3  area  1  rom  L.j  -uppe  r 
atld,  res  per  t  i  vr  I  V  , 

1  .or  (Trapezoid  2  area  from  I-'  .-lower) 
or 

Lot;  (Trapezoid  4  area  from  K.^. -upper) 

1  sine  tlie  data  sc  in  1  able  2  to  provide  the  desired  variates,  we  obtain*’. t  the 
■  oeffieient  s  given  in  f  able  5(b).  Considerations  of  I  able  5(a)  and  5(b)  show  that 
the  three -variate  regression  gives  a  better  lit  to  the  data  set  than  the  two-variat" 
regression.  The  higher  frequencies  also  are  better  it  titan  the  lower  irepucni  irs 
as  determined  t>v  the  correlation  coefficients,  to  the  >  30  IMeY  proton  data  than  to 
tlic  >  It)  Met  proton  data. 


1.5  tom  \  ;i null'  l  incuv  Ki’tiri'ssiun 

1  he  three  hit’ll  frequencies  of  llte  five  considered,  which  tend  to  five  better 
linear  regressions  to  the  proton  data,  were  broken  apart  again  'o  five  a  lour- 
'  .iriate  linear  regression  ol  the  lorm 


1  ’  1 


a...\.,  , 


Y  l  os  (ad;  list  ed  nr- itnn  1  ux ) 

>,  |  I  ,og  |  area  of  2(1. ':>  -MHz.  rectangle,  <  fi gure  '■  i j 

\.,  Lot;  Lire, i  of  t"l'5-Mllz  rectangle,  (figure  t'l| 

1 .  *  >  e  j  T’.oi)  -MHz  ,.,.et  int;le.  f  1  L  .  •  >  i  r-  '*r 

1  he  ire  is  >t  '  lie  t  luce  re  tanples  a  re  foun  1  hv  t  lie  n  "o  hr  '  ot  •  ’ ■■  into.,  rate  | 

radio  'lux  al  e.ieh  frcpienev  limes  ihe  wi  Itli,  in  hertz,  c  die  ••  •  ,'.T  i.  [hi  pro 

v  ides  a  measure  if  the  energv  in  a  frequenev  hand  c  mr  I  ■■  r  r  • ;  •  ■  • ; u« •  v  r  ih 
"nation.  I  rom  figure  h.  the  wid'h  of  the  rectangle-'  were  it,'  cm-  i  ■:  • ,  n  ,u 

I  lie  :*eoi "» -l  r  i  i  mean  '  i'''Mie'i  m-  w  .  :  -  »ls  om  i  hi  .  i  •> 

w  nirh  a  ;  ime  inlei;ra:e  i  'lu\  w.i,-.  •  miii  1  an  I  'die  i  i  \  e  "ir.  ' :  ciu-  ;;  ;•  ■  ,  I ,  ■  •  ,  a  |  t  „  1 . 

It  .  t  a  I  I  Ule.se  '  reoneiii  Il'S  !  lie  liolli'  I . ’ .  '  !'*  p '  |e  I!  ■  ill  i  'he  1  ill.'!  "  nil  :  •  I  Ueu 
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Figure  0.  Three; -Frequency  S|«  '  trum  of  Time- 
Integratod  Radio  Fluxes  of  L -portion  of  28  August 
1966  Solar  Radio  Burst.  Associated  rectangles 
used  in  four-variate  linear  regression  are  indicated 


The  difference  Between  the  upper-mean  frequency  and  the  lower-mean  frequency 
was  then  taken  as  the  width  of  the  frequency  band  around  each  observed  frequency. 
The  upper  frequency  used  to  determine  the  upper-mean  frequency  for  the  8800-MHz 
area  was  1.6,  100  MHz  which  is  the  next  highest  frequency  observed  at  the  Sagamore 
Hill  Radio  Observatory  and  the  other  observatories  of  the  Air  Weather  Service's 
Radio  Solar  Telescope  Network  flt.STN).  The  best  fit  linear  regression  coefficients 
were  then  obtained  using  the  data  in  fable  2.  Table  6  gives  the  results  of  this  re - 
L1  cession.  Tli"  correlation  coefficient,  R  ^  18  ‘  ailed  the  multiple  correlation 

coefficient  of  variate  '.  on  variates  2,  2,  and  -1.  These  correspond  to  V,  \ ^ ,  V,,  a 
res  iicrt  ive|  y,  m  the  expression  for  the  linear  lit  6'  a()  ■  a  X  j  •  <  j -V,  •  a^X.,.  I  he 

quantity  Sj  j  is  th'-  stand.tro  error  of  estimate  of  variate  1  on  variates  2,  i, 
and  1.-  '  As  would  be  collected,  the  correlation  coefficient  which  is  a  m<  astire  of 
'he  goodness  ol  the  'it  ol  ttic  cxnression  to  the  lata,  in  each  case  ts  the  itiuhest  one 
■  ilit  a  i  tie  I  t  n  t  h  i  s  -•  ud  v.  I  lie  •  .'Ul-Mc\  t  trot  on  data  t  re  al  so  seen  •  o  :  i  t  bet  t  e  r  t  ban 
'lie  •  It)  \|i-  i .  1 1  :  i  . 
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CONC.I.l  SIONS 

Sniiii'  considerations  have  hern  made  here  of  the  relationships  between  the 
integrated  cad'o  fluxes  of  solar  flares  and  the  associated  energetic  proton  flux 
increases.  Adjustments  to  the  particle  flux  for  the  position  :> f  the  flare  on  the  face 
>f  the  sun  have  been  applied  and  various  combinations  and  weightings  of  the  inte¬ 
grated  radio  fluxes  have  also  been  tried.  In  the  end,  if  tne  correlation  coefficient 
between  the  variatc-s  considered  is  a  measure  of  their  predictive  capability,  there 
appear  to  he  several  combinations  and/or  adjustments  which  are  of  comparable 
apabilitv.  Just  which  one  will  be  list'd  will  be  up  to  the  user  and  the  degree  of 
itiv olvement  desired. 

Predictions  of  impending  occurrences  attempt  to  foretell  what  will  take  place 
before  thev  actually  do  and  what  the  magnitude  of  the  occurrences  will  be.  The  re¬ 
sults  of  this  study  mav  he  of  use  in  the  following  situations: 

1.  I  he  influx  of  energetic  protons  into  the  polar  ionosphere  of  the  earth  causes 
an  increased  absorption  of  III  radio  signals  that  pass  through  the  polar  ionosphere. 

1  his  absorption  is  known  as  polar  cap  absorption  or  PCA.  I  he  magnitude  of  the 

absorption  can  be  related  to  tin-  it)  MoV  proton  flux  bv  the  following  expression 
2  0-1 
o  :  10. A  .  where  1’  is  the  proton  llttx  in  protons  (cm  sec  ster)  and  A  is  the 

O.J 

Jd-.MI1/  rtomi'tcr  absorption  in  I )  It.  ”  1  he  integrated  radio  llttx,  predicting  the 

:  •  i  a  xi  ftit' ,  proton  flux  increase,  will  make  it  possible  to  determine  the  absorption 
'  ( )  hi  •  i  •  \  p> '  ■  ted. 

2  1.  t  . t 1 1 ■  1 1 1 ,  .1.1’.,  and  Tarns!  rom,  (1.1..  (  1  9711 1  A  Catalog  of  Proton  Kver.ts 

I  oiik  -  I  'fi7 1  i|,o  jug  Non-Classical  Solar  Kadlo  Hurst  Spectra, 

\  IdTT  -TR  -  TTt  T)  12  t ,  A I )  A 0608  1  (T,  Knvi  ronmcntal  lies .  Pape  rs  622, 
llin-.oi:  \l  li,  .\l  assaehuset  t  s . 


2.  The  increased  incidence  of  energetic  protons  into  the  polar  ionosphere 
causes  the  height  of  the  D- region  to  be  lowered  bv  changing  the  D- region 
ionization  structure.  It  has  been  shown  that  during  a  proton  flux  increase,  the 
height  of  the  south  polar  D-region  tracks  the  variations  in  t he  proton  flux  reasonably 
well.-  Changes  in  the  l)-region  height  influence  the  propagation  of  VL.K  radio 
signals  across  the  polar  caps  causing  variations  in  their  phase  and  amplitude 

2  g 

structure,  which  impacts  their  use  for  radio  navigation  purposes.  The  integrated 
solar  radio  flux  should  make  it  possible  to  make  an  earlier  determination  of  the 
maximum  l)-region  height  change  to  be  expected,  from  which  the  disruption  of  the 
trans-polar  VLK  propagation  may  be  anticipated. 

3.  One  requisite  of  a  prediction  process  is  that  it  provide  the  prediction  as  far 

27 

ahead  of  its  occurrence  as  possible.  Smart  and  Shea,  suggested  a  method  of 
predicting  the  end  of  a  proton  flux  increase  by  extrapolating  the  end  time  of  the  in¬ 
crease  after  the  maximum  of  the  increase  has  occurred.  The  process  samples  'he 
proton  flux  levels  at  a  succession  of  times  to  obtain  the  time  character  of  the  flux 
drop-off  and  then  projects  the  time  when  the  flux  increase  will  be  reduced  to  zero. 
The  same  procedure  can  be  applied  to  the  radio  flux  increase  being  observed,  to 
predict  when  the  increase  returns  to  zero  or  to  some  predetermined  level.  With 
that  time  and  the  sense  of  the  functional  nature  of  the  radio  flux  decrease,  the 
time -integration  of  the  radio  fluxes  can  be  accomplished.  With  these  integrations, 
the  expected  peak  energetic  proton  flux  can  be  predicted.  This  prediction  mav  not 
be  the  best  possible  because  the  position  of  the  flare  on  the  face  of  the  sun  mav  not 
be  known,  but  it  will  at  least  give  a  first  approximation  of  the  proton  flux  increase 
to  be  expected.  por  a  radio  burst  of  comparatively  short  duration,  the  time  advan¬ 
tage  of  the  prediction  mav  not  be  great,  but  for  radio  events  with  longer  decav  times, 
and  for  long  duration  events  of  the  "nonclassical"  type  considered  bv  Castelli  and 
Tarnstrom”  the  timeliness  of  the  neak  proton  flux  prediction  mav  he  enhanced  bv 
as  much  as  an  hour  or  more. 


2o.  Helms,  W..I.  (HITS)  Polar  l)-region  electron  density  profiles  during  a  solar 
proton  event.  Radio  Sci.  I3(No.  5):8a3. 
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26.  Westerlund,  S.  ,  Reder,  1- .  11.  ,  andAbom.  UhfiO)  Kffeets  of  polar  cap  absorp¬ 

tion  events  on  YI.K  transmissions.  Planet.  Suttee  Sri.  l7(No.  7  h  132b. 

- — - : -  l/W 

27.  Smart,  I).  I- .  ,  and  Shea,  M.A.  (1h73)  Prediction  >f  the  end  of  solar  proton 

events.  Prof,  loth  Plenary  .Meeting  of  COS  PAR.  Madrid,  Spain, 

10-14  .Mav  P'72,  Space  Res.  \IIlT  \  <>1.  T,  Akadrm ie  ver lag.  Berlin. 
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